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>>■ ABSTRACT 

We use deep observations obtained with the Photodetector Array Camera & Spectrometer (PACS) onboard the Herschel space ob- 
servatory to study the far-infrared (FIR) properties of submillimeter and optically faint radio galaxies (SMGs and OFRGs). From 
C*^ ' literature we compiled a sample of 35 securely identified SMGs and nine OFRGs located in the GOODS-N and the A2218 fields. 

This sample is cross-matched with our PACS 100 /jm and 160 yum multi-wavelength catalogs based on sources-extraction using prior 
__ . detections at 24 yum. About half of the galaxies in our sample are detected in at least the PACS 160 yum bandpass. The dust temper- 

Q^ ' atures and the infrared luminosities of our galaxies are derived by fitting their PACS and SCUBA 850 yum (only the upper limits for 

r\ ' the OFRGs) flux densities with a single modified if} = 1.5) black body function. The median dust temperature of our SMG sample 

• . is Trfaj, = 36 ± 8 K while for our OFRG sample it is Tj^st = 47 ± 3 K. For both samples, median dust temperatures derived from 

r^^ . Herschel data agree well with previous estimates. In particular. Chapman et al. (2005) found a dust temperature of Trf,,.,, = 36 ± 7 K for 

Oh' a large sample of SMGs assuming the validity of the FlR/radio correlation (i.e., q = /ogio(Lf/K[W]/Li 4ghz[WHz"']/3.75 x 10'^)). 

The agreement between our studies confirms that the local FlR/radio correlation effectively holds at high redshift even though we 
find {q) = 2.11 ± 0.19, a slightly lower value than that observed in local systems. The median infrared luminosities of SMGs and 
OFRGs are 4.6x 10'^ Lq and 2.6x 10'^ Lq, respectively. We note that for both samples the infrared luminosity estimates from the radio 
^ ' part of the spectral energy distribution (SED) are accurate, while estimates from the mid-lR are considerably (~ x3) more uncertain. 

' Our observations confirm the remarkably high luminosities of SMGs and thus imply median star- formation rates of 960 M© yr"' for 

SMGs with 5(850/jm) > 5mJy and 460 Mgyr"' for SMGs with 5(850//m) > 2mJy, assuming a Chabrier IMF and no dominant 
AGN contribution to the far-infrared luminosity. 
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1. Introduction merging sytems (i Tacconi et"an l2006l l2008h with high star- 

forma tion efficiencies com pared to typical galaxies of a similar 



Herschel observations probe the rest-frame far-infrared emis- ^^^^ dPaddi et all llOOSh . Therefore, these SMGs are thought 
g sion of high-redshift galaxies. Thus, they provide for the first ^^ g^hibit very intense (SFR~1000 M^ yr') short-lived star- 
^ ;time robust estimates of the infrared luminosities of these formation bursts triggered by mergers and to be the high-redshift 
j^ . high-redshift galaxies and test previous measurements that were progenitors o f local m assive early-type galaxies (IDaddietall 
^ ^ based on extrapolation from shorter or longer wavelengths. We 2007a b: Tacconi et al.Ll200llci^atti et al.. 2008h. 
k^ here focus on two populations of high-redshift star-forming Although SMGs provide a powerful tool to constrain the 
^> ^ galaxies selected at submillimeter (submm) and radio wave- formation and the evolution of high-redshift dusty star-forming 
C^ I lengths. galaxies, their selection is subject to strong biases. In particular. 
Since their discovery in the late 1990s, submillimeter because submm observations probe the blackbody emission 
galaxies (SMGs) have become the selection of choice for the ^f dust in the Rayleigh-Jeans regime, they are strongly anti- 
most luminous tail of the high-redshift star-forming galaxy correlated with the dust temperature (S^so oc Tf'). For a given 
population It has been foundjhat SMGs have ajXEical redshift ^^f^^^^ luminosity, galaxies with hot dust mTght fall below 
of z ~ 2 ( iChapman et al-l |2005t Pope et aLj | 2006|), ar e mas- j^e detection limit of current submm instruments. The first 



P 'inrtA'- 



sive systems (M^ 10'" - 10" Mq, Swinbank et al., 2004^, observational evidence of a missing population of high-redshift 
Taccometal., |2006|) and are compact (e.g., Taccometal., dusty star-form i ng gal axies with hot dust has been given by 



200i. Interferometric observations of their CO molecular gas jchat^man et all dfool using a selection of radio-detected 



suggest that the most luminous SMGs {Ss50^n< > 5 mJy) are ^ut submm-faint galaxies with UV spectra consistent with 

high-redshift starbursts. These optically faint-radio galaxies 

Send offprint requests to: B.MagneWi, e-mail: Kagnellimpe.mpg.de (OFRGs) have a comoving volt ime d ensity (i.e., ~ lO^^Mpc"^ 

* Herschel is an ES A space observatory with science instruments pro- at 1 < z < 3, Chapman et al. l2004l), stellar masses and sizes 

vided by European-led Principal Investigator consortia and with impor- comparable to SMGs, and some have a dust temperature of 
tant participation from NASA. 
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Fig. 1. Spectral energy distribution (SED) of one SMG (left) and one OFRG (right). Red diamonds present our PACS measurements wliile green 
squares present the multi-wavelength ancillary data taken from the literature (Pope et al., 2006; Casev et al., 2009a b). The modified blackbody 
emission (yS = 1.5) best-fitting the data are shown by solid lines. Dashed lines present the Dale & Helou SED template best-fitting the mid- to 
far-infrared observations. In the left panel, the blue dotted-dashed line shows the Dale & Helou SED template best-fitting the submm and radio 
photometries. The inset in each panel shows ;^'- vs Tj„„. 



~ 52 K (ICasevetal.L l2009a''bS. 

While SMGs and OFRGs are an important component 
of the high-redshift massive galaxy population, many of their 
fundamental properties still rely on indirect measurements. 
In particular, direct determ i nation s of SMG dust temperatures 
were limited (iKovacs et al.L l2006l) because they ■were not done 
using rest-frame far-infrared observations on both sides of the 
peak of the SEDs. More importantly, their infrared luminosites 
are still debated. Indeed, because theoretical simulations of 
galaxy evolution have had great difficulties to account for the 
cuiTent inferred luminosities/star-formation rates and number 
counts (Baugh et al., 2005; Dave et a l., 2010), they still question 
whether these luminosities have been overestimated or whether 
the IMF is significantly more top-heavy than in the local 
Universe. 

Using deep observations by the Photodetector Array 
Camera & Spectrometer (PACS; Poglitsch et al. 2010) onboard 
the Herschel space observatory (Pilbratt et al. 1201 Oh we will 
have for the first time robust estimates of the dust temperatures 
and the infrared luminosities of SMGs and OFRGs. Throughout 
the paper we use a cosmology with Hq - 70kms"'Mpc"', 
Qa = 0.7 and Qm = 0.3. 



2. Observations 

In this study we used deep PACS 100 ;um and 160 /im 
observations of the Great Observatories Origins Deep Survey- 
North (GOODS-N; 12''36"', +62° 14') and the Abell 2218 
(16''35'", +66°12') fields. These observations were taken as part 
of the PACS Evolutionary Probe (PEFQ) guaranteed time key 
program. The GOODS-N field covers a region of 10' x 15' (30 
hours), while the deep part of the A2218 field covers a region of 
4' x4' (13 hours). 

At the resolution of Herschel, all sources in our fields are 
point sources (i.e. FWHM- 8" [12"] at 100 yum [160 ^im\). 
Flux densities are hence estimated using a point spread function- 
fitting technique based on prior source positions detected at 
24 yum. The use of priors provides a straightforward association 
between the IRAC, MIPS and PACS sources. Using Monte 



Carlo simulations we estimate the quality of our PACS 100 yum 
and 160 //m catalogs, i.e. photometric error, completeness 
and contamination as a function of the flux density. In the 
GOODS-N field our observations reach a 3cr limit of ~ 3 mJy 
and ~ 5.7 mJy at 100 yum and 160 /im respectively, while in the 
A2218 field they reach a 3cr limit of ~ 2.5 mJy and ~ 5 mJy at 
100 yiim and 160 yum respectively. 

A complete description of PEP data reducti on and sources 
extraction is given in Appendix A of Berta et al. (l2010h . 



3. Galaxy sample 

To obtain a robust measurement of the dust temperature and 
infrared luminosity of a given galaxy one needs to have an ac- 
curate estimate of its redshift. Consequently, we decided to re- 
strict our study to a sample of SMGs and OFRGs with accu- 
rate redshift estimates derived from secured radio/mid-infrared 
identifications (PACS identifications of SMGs are presented 
in Dannerbauer et al. in prep). In the GOODS-N field, our 
SMG sample is based on multi-wavel ength ident i ficatio ns of 
SCUBA and AzTE C sources made by iPope et all (l2006i PI and 
IChapin et al.l (l2009h . respectively. SMGs with tentative redshifts 
determined from their IRAC or mid/far-infrared/radio colors 
were excluded from our sample. Sources with multiple optical 
counterparts (GN04, GN07, GN19 and GN39) were treated as 
a single system (i.e. we will use the sum of the radio and mid- 
infrared flux from the two components when determining their 
far-infrared prope rties) because the y are all thought to be in- 
teracting galaxies dPope et alll2006h . All these different criteria 
yield an SMG sample containing 29 sources in the GOODS-N 
field. In the A2218 field, our SMG sample is assembled from the 
literature (Kneibetal., 2004; Knudsen et al., 2006, 2008) and 
contains six lensed sources. Because these galaxies are magni- 
fied, their mid-to-far infrared fluxes were de-magnified prior to 
further analysis using magnification factors from the above refer- 
ences. Among these six lensed sources, three coiTespond to the 
same lensed galaxy (SMMJ16359+661 2; Kneibetal., 2004 . 
Finally, our OFRG sample is taken from lCasev et al.l (l2009allbh 



http://www.mpe.mpg.de/ir/Research/PEP 



- For GN05, GN07, G NIO, GN20 and GN 20. 2, we used the spect ro- 
scopic redshifts revised in lPope et alJd2008l) and lDaddi et alJ J2009al lbh. 
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Fig. 2. Dust temperature-luminosity relation. The filled blue squares 
and opened blue squares denote SMGs located in the GOODS-N fields 
with spectroscopic and photometric redshift, respectively. The crosses 
denote sources which contain an AGN, as indicated by the presence of 
hard X-rays. The filled green diamonds present SMGs located in the 
A2218 field (the three data points with Lm ~ 1 x 10'' L© correspond to 
the same lensed galaxy). The filled red circles are from our OFRGs 
sample. The striped area presents results for SMGs extrapolated by 
[Chapman et al. ( 2005) from radio and submm data. The lChapman et al] 
(2003) derivation of the median and interquartile range of the Tjust-Lm 
relation observed at z ~ is shown by solid and dashed-dotted lines, 
linearly extrapolated to 10'^ Lq. 



and contains nine sources, all situated in the GOODS-N field. 
We note that all but three sources of our entire sample (i.e. SMGs 
and OFRGs) have spectroscopic redshifts. 

The SMG and OFRG samples were cross-matched with our 
R4CS multi-wavelength catalogs using a matching radius of 3". 
We detected 19 out of 35 SMGs in at least the PACS 160 yum 
bandpass (17 out of 33 if not multi-counting the 3-component 
lensed source detected in A2218). The PACS sample is slightly 
biased towards lower redshift sources because of the positive K- 
correction : while the median redshift of our parent SMG sample 
is z = 2, the median redshift of our PACS detected SMG sam- 
ple is z = 1 .7. Five out of nine OFRGs have PACS 100 jim and 
160 jiva detections. This sample is also slightly biased toward 
lower redshift {z - \ .5) because the OFRG situated at the high- 
est redshift is undetected in our PACS images. 

We note that our SMG sample contains sources with 
2 mJy < S 850;/m < 5 mJy, while the most luminous tail of the 
SMG population, mostly associated with major mergers, is de- 
fined using S'sso^m > 5mJy. Below, we will draw our conclu- 
sions distinguishing these two populations of SMGs. 



4. Data analysis 

In order to infer the dust temperature of our galaxies we fitted 
their PACS and SCUBA photometry (only the upper limit for 
the OFRGs) with a modified blackbody function, with a dust 
emissivity jS =1.5 (see Fig[TJ. Their total infrared luminosities 
(L/R[8-1000|/m]) were infeiTed from these best fits using the 
far-infrared lum inosit y definition (L/r/R[40 - 120/ im]) given 
by He lou et al. (Il988h and a color-correction term dPale et al.L 
|200ll^/fi = l-91xL„R). 

We adopted a single dust temperature characterization 
because studies of IRAS galaxies have demonstrated that this 
provides an accurate diagno stic of the ty pical heating condition 
in their interstellar medium ( iDesert et al.L 11990). While for most 
of our galaxies this single dust temperature characterization 
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Fig. 3. Submm flux densities as function of the infrared luminosity. 
The symbols are same as in Fig|2]except for the OFRGs, for which we 
have only upper limits. The solid and dashed lines show the linear fit 
to the Ss5o^m-LiR relation and the Icr envelop (L;k[Lo] = lo"-^-"^ x 
^Isopm [mJy])- If we remove the lensed-SMGs of A2218 from the fit, 
we find a weaker correlation (L,h[Lq] = W^oo+oa ^ 5i5±ai [mjy]). 



provides a good description of their far-infrared SED, for 6 
SMGs, this single dust temperature model yields high;^'^ values 
(i-S-X^ > 2.7 1 + Ndof)- All these galaxies appear either to be the 
more distant ones or to exhibit far-infrared colors typical of very 
cold systems. In both cases, their PACS 100 /im flux densities 
might by contaminated by a hotter dust component. Therefore 
we excluded their PACS 100 //m photometry from the fit and 
recomputed their dust temperatures. We note that excluding 
the PACS 100 i-im photometry from the fit of all our galaxies 

changes their median dust temperature by only ATjusi 2 K. 

We tested that our results are relatively insensitive to j3 
and to our single dust component characterization. Indeed we 
note that using /? = 2, we found only small differences in the 
values of Tdusi (ATdust ~ +3 K). Moreover, we also note that 
to fit the PACS and SCUBA photomet ry with the multiple dust 
components model of Dale & Helou (l2002h yields rest-frame 
infrared colors (SeOfjm/S \oo^lm), or equivalently Tdusi, which 
excellently agree with those inferred from our blackbody 
analysis. 



5. Discussion 

Figure |2] depicts the locations of our SMGs and OFRGs on 
the Tdiisi-LiR plane and Fig |3] shows their locations on the 
S&sOfim-LiR plane. As already mentioned, OFRGs are biased 
towards hot dust temperatures; their median Tdusi is of 47 + 3 K 
and their median L/r is 2.6 x IO'^Lq. In contrast, SMGs have 
lower dust temperatures with median Td„si - 36 + 8 K and 
LjR = 4.6 X 10'^ Lq. We note that lensed-SMGs from A2218 and 
with L/R < 2 X 10'^ Lq exhibit intermediate dust properties and 
are less biased towards cold dust temperatures than the entire 
SMG sample. This is because these galaxies would have escaped 
both the SMG and OFRG selection method without magnifi- 
cation. We also note that bright SMGs (i.e. ^gso/jm > 5mJy) 
have higher median infrared luminosities (L/r = 9.6 x 10'^ L©) 
and higher median dust temperatures (Tdust = 38 K) than the 
entire SMG sample because there is a correlation between 
Sssoum and Ljr (Fig [3]l. These estimates are the first direct 
observational measurements of the dust temperatures and the 
infrared luminosities of SMGs and OFRGs. 

Our observations reveal that high redshift dusty star-forming 



Magnelli et al.: Far-infrared properties of SMGs and OFRGs 



galaxies exhibit a wide range of dust temperatures. In particular 
at low infrared luminosities (L/r < 4 x IO^'Lq) the dust temper- 
ature dispersion observed in our sample might suggest a higher 
Tdust-LiR scatter than that observed by IChapman et alJ (l2003h 
at z ~ 0. Nevertheless, this conclusion is most likely driven by 
selection effects because a significant fraction of the galaxies 
with intermediate dust properties were probably missed by our 
current sample. Indeed, we note that studying a L/^-selected 
sample of galaxies observed with Herschel, Hwang et al. (in 
prep) find modest changes in the Tdust-Lm relation as function 
of the redshift : at z > 0.5, galaxies with L/r > 5 x IO'^Lq are 
slightly colder (~ 3 K) than local ones and the scatter of the 
Tdiist-LiR relation slightly increase at high redshift. 

Though previous estimates of the dust temperatures 
of SMGs and OFRGs reUed on indirect observations, they 
agree relatively we ll with our measurements. In partic- 
ular, IChapman et al.l (l2005h found a dust temperature of 
Tdust = 36 ± 7 K for a large sample of SMGs assuming the 
vaUdity of the FIR/radio correlation. In order to establish this 
agreement on a common sample, we applied the same method 
as Chapman et al. (2005) to our SMG sample, i.e. we fitted 
the radi o and 850 um pho tometries with dust SED templates 
from .Dale & Heloul (l2002h and then translated them into Tjus, 
using their R(60,100) to Tjust map. With this method we found 
higher Tdust (by ~ 4K) and Lju (~ xl.5 times) than what we 
obtained using our blackbody analysis. These discrepancies 
arise because the Dale & Helou SE D templa t es ass ume a 
FIR/radio correlation witii (q) = 2.34 dYunetali 1200 ih . while 
in ou r samples we find (q) = 2.17 + 0.19 (see also Ivison et al. 
1201 Q l, this issue). Although this value of (q) is still in line with 
results from local systems (which have a dispersion of 0. 19 dex), 
it is also consistent with an evolution of (17) proportional to 
(1 +2)-oi5±o.o3 ^^ found bv ilvison et all dloia) . 

Our results reveal that one can obtain a very reliable 
estimate of the infrared luminosity of a given galaxy from 
its radio flux density {a-[logm{L'^f" llf^''"'"'''')] -0.18 dex) 
using {q) - 2.17. In contrast, the use of the 24 yum emission 
and of the Chary & Elbaz SED library yields an inaccurate 
estimate of the infrared luminosity characterized by a large 
scatter (cr[logio(Ljj^'^'" /Lj^'^ " ■*)] ~ 0.48 dex) and a systematic 
overestimation (~ x2 times) of the most luminous galaxies 
(LiR > 4 X lO'^Lo). We thus find that for the very luminous 
infrared galaxies studied here, luminosity extrapolations based 
on the radio emission are considerably more reliable than those 
based on the mid -infra red emission (see also Elbaz et al. 1201 Ol 
and Nordon et al. l2010l) . 

Using (q) - 2.17, we can predict the infrared luminosities 
of our PACS undetected SMGs. Then, using the Dale & Helou 
SED templates normalized to these infrared luminosities, we 
can fit their 850 //m photometries. We find that for 13 out of 16 
undetected SMGs, PACS flux densities inferred using these fits 
are below the detection threshold of our observations. Of the 
three sources with PACS fluxes predictions above our detection 
threshold, one is known to be contaminated by an AGN and 
the other two are suspected to have wrong redshift estimates 
(iDaddi et all '2009b). The PACS nondetections are thus fully 
consistent with the properties inferred from the detections. This 
analysis cannot be performed on our PACS undetected OFRGs 
because they are also undetected at the SCUBA wavelength. 

Our observations unambiguously confirm the re- 
markably large infrared luminosities of bright SMGs (i.e. 
Ss50^m > 5mJy) which correspond to SFRs of 960 Mq yr ' 



(SFR [Mo yr-i] = 1 x IO^^^Lir [Lq], assuming a Chabrier 
IMF and no dominant AGN contribution to the far-infrared 
luminosity). Such hig h SFRs are diffi cult to reconcile with 
secular evolution (e.g. iDave et all l2010h and could correspond 
to a merger-driven stage in the evolution of these galaxies. This 
hypothesis is supported by CO observations of bright SMGs 
which have re vealed large C O line-widths and disturbed gas 
morphologies dTacconi et al.L |T008). Our observations also 
confirm that OFRGs exhibit higher dust temperatures than 
faint SMGs (8^50 f^^ < 5 mjy) observed at the same redshift 
and with equivalent infrared luminosities. While the relatively 
low median SFRs of OFRGs and faint SMGs (260 Mq yr ' 
and 106 M© yr"' respectively) could be explained by secular 
evolution, we need to understand why they exhibit different dust 
temperatures and to study a possible link with bright SMGs. A 
clear evolutionary picture will require detailed studies of the 
dust and molecular gas distribution in a sample of high-redshift 
star-forming galaxies unbiased towards any particular dust 
temperature. This SFR-selected sample can now be built with 
our ongoing deep Herschel observations. 
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